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ABSTRACT 

In the mammalian mitochondrial translation appar- 
atus, the proteins and their partner RNAs are coded 
by two genomes. The proteins are nuclear-encoded 
and resemble their homologs, whereas the RNAs 
coming from the rapidly evolving mitochondrial 
genome have lost critical structural information. 
This raises the question of molecular adaptation of 
these proteins to their peculiar partner RNAs. The 
crystal structure of the homodimeric bacterial-type 
human mitochondrial aspartyl-tRNA synthetase 
(DRS) confirmed a 3D architecture close to that of 
Escherichia coli DRS. However, the mitochondrial 
enzyme distinguishes by an enlarged catalytic 
groove, a more electropositive surface potential 
and an alternate interaction network at the 
subunits interface. It also presented a thermal sta- 
bility reduced by as much as 12 C. Isothermal titra- 
tion calorimetry analyses revealed that the affinity 
of the mitochondrial enzyme for cognate and non- 
cognate tRNAs is one order of magnitude higher, 
but with different enthalpy and entropy contribu- 
tions. They further indicated that both enzymes 
bind an adenylate analog by a cooperative allosteric 
mechanism with different thermodynamic contribu- 
tions. The larger flexibility of the mitochondrial 
synthetase with respect to the bacterial enzyme, in 
combination with a preserved architecture, may rep- 
resent an evolutionary process, allowing nuclear- 
encoded proteins to cooperate with degenerated 
organelle RNAs. 



INTRODUCTION 

Mitochondria are of endosymbiotic origin and have 
undergone massive gene transfer to the host nucleus 
during evolution (1). On sharp reduction (~70-fold) to a 
16.6 kb circular mammalian mitochondrial (mt) DNA, 
only 37 genes remain coding for 13 respiratory chain 
subunits, 22 tRNAs and 2 rRNAs (2). The faster evolu- 
tionary rate of mt-DNA than the nuclear genome (3,4) 
leads to abnormal RNAs, shrunken in size and often 
lacking important signals. mRNAs are deprived of 
3'- and 5'-untranslated regions, rRNAs are smaller than 
their bacterial or eukaryotic homologs, and an RNA 
moiety for RNase P is absent (5-7). Most mt-tRNAs 
have shortened sequences, miss crucial folding and recog- 
nition nucleotides as compared with 'classical' tRNAs and 
are more flexible (6,8-10). Hence, preservation of mito- 
chondrial translation relies on a continued adaptation of 
proteins to the peculiarities of partner RNAs. One com- 
pensatory mechanism consists in an increase in the 
number and/or size of partner proteins (encoded in the 
nucleus and imported) that take over the role of missing 
RNA domains. For instance, bovine mt-ribosomes 
contain 30 proteins in the small subunits and 48 in the 
large subunits, instead of 21 and 33, respectively, in pro- 
karyotes, with N- or C-terminal extensions (11,12). 
Human RNase P is fully proteinaceous with three 
protein subunits (7). 

Mechanisms to compensate for degeneration of 
mammalian mt-tRNAs remain mainly unsolved, espe- 
cially for mt aminoacyl-tRNA synthetases (mt-aaRSs; 
specificity indicated by the amino acid abbreviated in a 
one-letter code, the origin of the aaRS is indicated by a 
three-letter abbreviation of the organism's name with the 
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mitochondrial origin given by the digit 2 e.g. HsaDRSl 
for Homo sapiens mitochondrial aspartyl-tRNA 
synthetase), the enzymes which activate and bind a 
specific amino acid to their cognate tRNA (13). These 
enzymes are nuclear encoded and, with only few excep- 
tions, are of about the same size, have the same structural 
2D organization and contain the expected specificity 
signals and synthetase class signatures as bacterial and 
cytosolic aaRSs (14-17). The first 3D structures of mam- 
malian mt-aaRSs confirmed the absence of additional 
domains (18-21). 

Properties of mammalian mt-aaRSs enabling them to 
deal with mt-tRNAs remain mainly unknown (9,22,23). 
On one side, their catalytic efficiency is usually one to 
two orders of magnitude below that of bacterial aaRSs, 
aminoacylation identity rules rely either on minimalist sets 
of determinants in the tRNA or even are disobeyed 
(22,23). On the other hand, a long standing remarkable 
feature of mt-aaRSs is their high capacity in tRNA 
accommodation: mt-aaRSs charge both mt- and 
bacterial tRNAs, whereas bacterial aaRSs cannot cross- 
aminoacylate mt-tRNAs (24). Non-aminoacylation of an 
mt-tRNA by a bacterial enzyme is linked to the absence 
of essential identity signals and to the large flexibility of 
the mt-tRNA (25). These data suggest the existence of 
peculiar dynamic properties for the mt-aaRSs, which are 
not shared with a bacterial homolog. 

Here, we compared several biophysical properties 
of human mitochondrial aspartyl-tRNA synthetase, 
HsaDRS2, with them to those of a bacterial (E. coli) 
homolog, £coDRS. The crystal structure of HsaDRS2 
was solved, thermal stabilities were quantified, and 
thermodynamic parameters of substrate binding were 
determined by isothermal titration calorimetry. As an 
outcome, both proteins share a common 3D architecture 
but the mammalian mt-aaRS gained in plasticity, enabling 
it to handle any tRNA Asp and especially the non- 
canonical mt-tRNA Asp . We propose that enlarged plasti- 
city represents a co-evolution mechanism of the nuclear 
and mt genomes in eukaryotic cells subjected to highly 
different evolutionary pressures. This plasticity will be 
an important parameter to consider for the molecular 
understanding of a rapidly increasing number of human 
disorders linked to mutations in mt-aaRSs. 



MATERIALS AND METHODS 

Protein and tRNA purification 

iicoDRS was overproduced in E. coli JM83, purified 
to homogeneity and concentrated to 40 mg ml -1 (25). 
HsaT>RS2 [mt-AspRS(I41) in (26); Supplementary 
Figure SI] was cloned into pDEST14 using the Gateway 
technology (Invitrogen) and overproduced in E. coli BL21 
(DE3) Rosetta 2 strain as described (26). It was purified by 
affinity chromatography on an Ni-NTA column followed 
by gel filtration and concentrated to 20 mg ml -1 (concen- 
tration of monomers) in 50 mM of HEPES-Na pH 7.5, 
150mM of NaCl, 1 mM of DTT, 0.1 mM of ethylenedia- 
minetetraacetic acid and 10% (v/v) of glycerol. This 
buffer, called herein 'isothermal titration calorimetry 



(ITC) buffer', prevents denaturation of the mt-aaRS for 
several days as verified by dynamic light scattering (DLS) 
(Supplementary Figure S2). Magnesium had to be kept 
out, as even low concentrations led to precipitation 
of HsaDRS2. It was also used for the preparation of 
EcoDRS and of the tRNAs. tRNAs (sequences are 
indicated in Supplementary Figure SI) were produced by 
in vitro transcription and were purified to single nucleotide 
resolution on denaturing 12% polyacrylamide gel electro- 
phoresis. They were solubilized in ITC buffer and 
denatured/renatured before use as described (16). 

Aminoacylation 

Assays were performed at 25° C on transcripts as described 
previously (16,27). Specific activities were measured by 
incubating 0.2-3 nM £coDRS or 77-930 nM #raDRS2, 
with an excess of 4 uM cognate or non-cognate transcript. 
Kinetic parameters k cat and K M were derived from assays 
containing 2 nM of £coDRS or 50 nM of HsaDRS2 with 
0.2^t uM of Eco tRNA Asp or mt-tRNA Asp transcript. 

Crystallography 

Crystallization conditions were searched at 20°C using 
commercial kits as described (26). HsaDRS2 was 
crystallized by vapor diffusion with 10% (m/v) 
PEG-3350 and 0.5 M of ammonium sulfate as crystallants 
in 25 mM of Bis-Tris pH 5.5 (26). Resulting needle-like 
crystals (~300 x 20 x 20 um 3 ) were cryocooled in their 
mother liquor. They were extremely radiation sensitive 
and diffracted weakly. Complete data sets were collected 
with a fast PILATUS 6 M pixel detector (DECTRIS) at 
the X06DA beamline, Swiss Light Source (Villigen). 
Diffraction data at 3. 7 A resolution in P2 space group 
were processed using XDS and XSCALE (28) and were 
converted into structure factors without any sigma cut-off 
using TRUNCATE from the CCP4 package (29). 
Molecular replacement solutions were sought with data 
truncated at 4.5 A using a homology model derived from 
the structure of £coDRS (30) and the PHASER program 
(31) as implemented in the PHENIX package (32). 
Two HsaDRS2 dimers were identified in the monoclinic 
asymmetric unit. Refinement using a dynamic elastic 
network (DEN) as implemented in CNS 1.3 (33) led to 
a reorientation of the insertions and to a dramatic im- 
provement of the electron density map. All residues of 
HsaDRS2 were built except the 26 C-terminal amino 
acids (including the 6-His tag), which are disordered. 
The final coordinates (residues 42-630 for each of the 
four monomers) and the experimental data have been 
deposited with the PDB (ID: 4AH6). Figures were 
prepared with Pymol (www.pymol.org), electrostatic 
surface potentials were calculated using its APBS plug-in 
and are depicted at the same color scale in all figures. 

Dynamic light scattering 

Particle size analysis was performed at 25°C using a 
Malvern Zetasizer NanoS instrument and a 40 ul quartz 
cell containing 20 ul of protein solution at 0.5-8 mg ml -1 
in ITC buffer. Five consecutive measurements were col- 
lected in automatic mode, and data were processed with 
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the manufacturer's software. The hydrodynamic diameter 
dh was derived from diffusion coefficient via the Stokes- 
Einstein relation and was corrected for solvent viscosity 
and refractive index. Thermal stability was estimated by 
subjecting protein at 3mg ml -1 to temperature gradients 
from 20°C to 80°C. The temperature at which particle size 
was increased by a factor 2 was taken as the melting point 
Tm. Additives were prepared in water. 

Differential scanning fluorimetry 

Samples made of 20 ul protein solution at 3mg ml -1 and 
0.5 ul of SYPRO Orange™ 5000 x stock diluted 20-fold 
in dimethlsulfoxyde (DMSO) were transferred in 0.2 ml 
tubes. Additives prepared as described earlier in the text 
were brought in 1 ul samples, so that their final concentra- 
tion was as indicated. The variation of fluorescence inten- 
sity during a temperature gradient from 25°C to 95°C was 
monitored in a Stratagene Mx3005P quantitative polymer- 
ase chain reaction instrument. On unfolding of the 
protein, hydrophobic residues become accessible to the 
dye, which becomes fluorescent (34). Melting point Tm 
was deduced from maximum of the derivative plot. They 
are ~5°C higher as compared with data obtained by the 
DLS method. 

Synchrotron radiation circular dichroism 

The experimental set-up at the DISCO beamline, SOLEIL 
synchrotron (Saint-Aubin, France) was calibrated for 
magnitude and polarization with a 6.1mg/ml 
D-10-camphorsulfonic acid solution. Ten micrograms 
per milliliter of DRS solutions in 100 mM of potassium 
phosphate, 50 mM of KC1, 10% (v/v) of glycerol and 
1 mM of Tris(2-carboxyethyl)phosphine hydrochloride, 
TCEP, were placed in a CaF 2 cuvette, with an optical 
path of 8 um. Three spectra between 1 70 and 280 nm 
were recorded at temperatures ranging from 20° C to 
80° C by 4°C steps to assess the thermal stability of 
EcoDRS and HsaDRS2. Data were processed (spectrum 
averaging and solvent baseline subtraction) using CDtool 
(35). 

Isothermal titration calorimetry 

Before ITC measurements, conditions were determined in 
which the four investigated macromolecules (HsaDRS2, 
EcoDRS, Hsa mt-tRNA Asp and Eco tRNA Asp ) were 
stable during the time of the experiment. This was found 
to be the case in the so-called 'ITC buffer' (50 mM of 
HEPES-Na pH 7.5, 150mM of NaCl, ImM of DTT, 
0.1 mM of ethylenediaminetetraacetic acid and 10% of 
glycerol). Magnesium had to be omitted, as it leads to 
precipitation of HsaDRS. To eliminate possible errors to 
ITC experiments, all samples (proteins and tRNAs) were 
extensively and extemporaneously dialyzed against a same 
stock solution of 'ITC buffer', degassed with argon and 
ultracentrifuged at 4°C. ITC experiments were performed 
at 25°C in a MicroCal iTC 2 oo instrument (GE Healthcare). 
In a typical experiment, 2 ul samples of 300 uM £coDRS 
or of 100 uM HsaDRSl were injected into 200 ul of 36 uM 
or 1 1 uM, respectively, cognate tRNA transcript solution. 
Injections lasted 4 s at 120 s intervals. Cross-bindings were 



also performed with protein/RNA concentration in a 
~10:1 molar ratio. Retrieved AH and AG values were 
used to calculate the binding entropy AS (expressed as 
-TAS, with AG = AH — TAS). The adenylate analog 
5'-(3-[jV-(L-aspartyl)sulfamoyl] adenosine (Asp-AMS) 
was titrated against 200 ul of enzyme with 0.7 ul injections 
to HsaDRSl at 32 uM and 1.5 ul to .EcoDRS at 45 uM. 
Data were analyzed with the Origin software. 



RESULTS 

Solution properties of HsaDRSl 

HsaDRS2 and IscoDRS are homologous dimeric enzymes, 
both of bacterial evolutionary type, sharing a same size 
and 43% sequence identity, including all typical signature 
motifs of DRSs and the same secondary structure organ- 
ization (16,36) (Supplementary Figure SI). Recombinant 
proteins were overexpressed in E. coli and the final enzym- 
atic preparations brought into the same buffer ('ITC 
buffer', see 'Materials and Methods' section) suitable for 
ultimate comparative isothermal titration calorimetry ex- 
periments. These conditions were a compromise for 
non-aggregation of HsaDRSl, the more sensitive of the 
two enzymes (27), preservation of activity for both 
enzymes and stability of the tRNAs. HsaDRSl studied 
herein is deprived of the 40 first amino acids likely encom- 
passing the mitochondrial targeting sequence (26). Both 
proteins had the same mean hydrodynamic diameter (d h ) 
of 10 ± 1 nm in DLS, in agreement with a homodimeric 
structure (Supplementary Figure S2). Kinetic parameters 
and catalytic efficiencies of both enzymes tested under 
identical conditions confirmed the ~60-fold higher 
aminoacylation rate and efficiency in the EcoDRS 
cognate system compared with HsaDRSl (Table 1). 
However, £coDRS was unable to aminoacylate the 
mt-tRNA, whereas the mt-aaRS charged both cognate 
and non-cognate tRNAs with similar efficiency. 

Crystal structure of HsaDRSl 

Single needle-like crystals of HsaDRSl were obtained in 
the presence of PEG-3550 and ammonium sulfate (26). 
They yielded a complete diffraction data set (26). In 
spite of their modest resolution (3.7 A), these data were 
suitable to determine the architecture of HsaDRSl (26). 
The structure, solved herein by molecular replacement 
using a homology model derived from the structure of 
the iscoDRS (30), confirmed the dimeric status of the 
enzyme. It clearly showed the anticodon binding and cata- 
lytic domains of the monomers but did not enable visual- 
ization of the bacterial insertion domain (residues 
337^170, Supplementary Figure SI). This suggested 
either a structural rearrangement or a local flexibility of 
this domain. The ambiguity was solved by using a 
DEN-assisted refinement procedure (33) that allowed re- 
positioning of the insertion in well-defined electron density 
(Supplementary Figure S3 and Table 2). This domain re- 
orientation, different from EcoDRS, was not an effect of 
crystal packing because the four independent monomers 
in the asymetric unit adopted the same conformation. 
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Table 1. Comparative kinetic properties of HsaDRS2 and &oDRS in tRNA aspartylation 



DRS a tRNA Asp (transcript) k cat (lO^sec" 1 ) K M (uM) k CM /K M (lO^sec" 1 uM _1 ) L AAG b (kcal moP 1 ) 

HsaDRS2 mt 4.0 ± 0.4 1.1 ± 0.2 3.7 ± 0.7 72 -2.45 

//.soDRS2 Eco 10.3 ±2.0 1.2 ± 0.1 8.3 ± 0.8 35 -2.04 

EcoDRS Eco 230 ± 6.0 0.8 ± 0.1 293 ± 14.4 1 

EcoDRS mt nd nd nd >10 6 

All data were obtained at 25°C (T = 298°K) and correspond to mean values from at least duplicated experiments with independent and freshly 
prepared enzyme and tRNA lots. 

a /c ca i and K M values were derived from assays containing 2 nM of &oDRS or 50 nM of HsaDRS2 and 0.2^1 uM of Eco or Hsa mt-tRNA transcript. 
L stands for the relative loss in aminoacylation efficiency referred to the tRNA aminoacylation efficiency of £roDRS for Eco tRNA Asp [L = (k cal / 
Km)&o system/(k ca i/KM)other]; bAAG is the variation of the free energy change at the transition state when comparing aminoacylation in the 
homologous Eco system with the reactions in the heterologous or Hsa mt-systems [AAG = — RT InL] (37); nd for not detectable. 



Table 2. Crystallographic analysis of HsaHRS2 



Data collection statistics 
Beamline 

Space group (number) 
Unit cell lengths a, b, c (A) 
Unit cell angles oc, p, y (°) 
Resolution range (A) 
Highest resolution shell (A) 
Number of observationsO 
Number of unique reflections 
Completeness (%) 
Multiplicity 

Rmerge/Rmeiis ( '°) 

<Ila(Ij> 

Crystal mosaicity (°) 
Asymmetric unit content 
Solvent content (%) 
Refined atomic structure 
Resolution range (A) 
R-factor/R-free (%) 
Number of protein residues/atoms 
Rmsd and bonds (A) and angles (°)° 
Ramachandran plot b residues in most 

favored unfavored regions (%) 
PDBid 



SLS/X10A 
P2 (3) 

142.4, 82.6, 146.3 
90, 100.4, 90 
82-3.7 
3.8-3.7 

255 882 (9662) 
35 768 (2403) 
98.8 (91.2) 
7.2 (4.0) 

18.0 (63.9)/19.4 (72.5) 

9.4 (2.1) 

0.25 

2 dimers 
62 

30-3.7 

22.1/28.0 

4x 589/18 840 

0.005/1.08 

89.2/0.2 

4AH6 



Values in parentheses are for the highest resolution shell. 
"Rmerge = E AW E, \h(hkl) - <1 ( hkl ) >|/E M/ E,- I,(hkl) and redundancy- 
independent; R m „ as = E /at/ (n/n-l) 1/2 E, \I,(hkl) - <I(hkl)>\/Y: hkl E, 
Ii(hkl). 

b Analysis carried out with Molprobity in Phenix package. 



Altogether, 96% of the sequence of the crystallized 
homodimeric HsaDRSl is seen (Figure 1A), with a 
slight asymmetry between its two subunits because of vari- 
ations in the bacterial insertion domain (Supplementary 
Figure S4 and Supplementary Table S5). Overall, the 
architecture of HsaDRSl is similar to that of bacterial 
DRSs, in agreement with sequence conservation (16,36). 
Comparison with £coDRS crystal structure (30) high- 
lighted a strong conservation. Respective DRS cores, 
including the anticodon-binding domain, which adopts 
the so-called OB fold (for oligonucleotide-binding fold), 
the hinge region and the catalytic domain with its charac- 
teristic antiparallel (3-sheet (30,38) displayed an over- 
all rmsd (root mean square distance) of ~2A 
(Supplementary Table S5). These three domains shared a 
similar position with respect to the dimer interface. The 
main difference between the two aaRSs resulted from the 



bacterial insertion that underwent a rigid-body counter- 
clockwise rotation of ~26° and clearly made the mito- 
chondrial active site more open for the binding of tRNA 
(Figure IB and C). The overall fold of the insertion 
remained unchanged. This reorientation explains why 
the insertion was barely visible on the initial maps. 
Minor conformational diversity occurred in protein 
segments that diverge in sequence, located in the periphery 
of the catalytic domain, the flipping loop and the motif 
2 loop (short insertions/deletions, see Supplementary 
Figure SI), which are extremely mobile in DRSs (39,40). 

Based on the structural knowledge on human mt- 
tRNA^ 15 (41) a model of the complex with HsaDRSl 
could be built (Figure 2) and allowed for comparison 
with the bacterial complex. Notable alterations in the 
charge distribution within the tRNA interaction surface 
of the proteins were observed. Electropositive patches 
render HsaDRSl more attractive for its RNA partner 
(Figure 2), similar to what was observed in the three 
other mt-aaRS crystal structures solved so far (18-21) 
(Figure 5). Ten additional basic residues (nine Lys and 
an Arg highlighted in blue in Supplementary Figure SI) 
are inserted in the close neighborhood of HsaDRSl cata- 
lytic cleft. They may contribute to make the surface more 
attractive or directly contact the tRNA backbone, espe- 
cially in the hinge region where such kind of aspecific 
interactions with the D-stem was described in the E. coli 
complex (30). Otherwise, almost all residues making direct 
and specific interactions with bases at both ends of the 
Eco tRNA are conserved in the human enzyme (indicated 
in green in Supplementary Figure SI). This is true for 
E270/R271/R274 or R599, which may contact C74 or nu- 
cleotides 69-70, respectively, as well as for Q129/E140, 
R75/Q91/R125, R125, which are adequately positioned 
to interact with the anticodon bases C34, U35 and G36, 
respectively. A few residues (indicated in orange in 
Supplementary Figure SI) differ without major implica- 
tion for the conservation of the interaction network. 
Thus, the main modification is the presence of a Gly at 
position 269, instead of the usually conserved Asp, which 
accounts for the insensitivity of HsaDRSl to the nature 
of the discriminator base 73 (42). 

The dimeric interfaces of both DRSs differ notably by 
the distribution and nature of subunit/subunit contacts 
(Supplementary Figure S6). In total, 70 hydrogen bonds 
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Figure 1. Newly solved 3D structure of HsaDRSl in comparison with 
£roDRS. (A) Cartoon representation of HsaDRSl dimer. Characteristic 
DRS domains are highlighted in the left monomer: anticodon binding 
domain in yellow, hinge region in orange, catalytic domain in red and 
bacterial insertion in violet. The right monomer is displayed in black with 
position of residues conserved in HsaDRSl, and &oDRS is colored in 
green. Sequence identity is 40% overall, 47% in the catalytic domain, 
41% for the anticodon-binding domain and drops to 25% in the inser- 
tion. (B) Superposition of DRS monomers: HsaDRS2 chain colored as 
in (A), E. coli free (PDBid: 1EQR) and tRNA-bound (PDBid: 1C0A) 
monomers in light and dark blue, respectively. The bacterial insertion 
in HsaDRSl is rotated by 26° compared with its position in £foDRS 
(see arrow). Bottom left: zoom on superposed insertions illustrating 
the overall fold conservation. (C) Zoom on the active site groove with 
same color code for DRSs as in (B), and Eco tRNA is represented 
as a semi-transparent orange backbone. The distance between two 
conserved residues (in green) on each side of the groove is displayed: 
Phe390 in HsaDRSl (Phe340 in &oDRS) in the external helix of the 



and 28 salt bridges stabilize the iicoDRS interface, 
whereas HsaDRS2 only contains 60 hydrogen bonds 
and 20 salt bridges. Considering a 10% larger interface 
area (5500 versus 5100 A 2 ), the contribution of specific 
interactions is 25% lower per A 2 , whereas the contribution 
of aspecific van der Waals interactions is increased, 
which may on average lead to a less cohesive HsaDRS2 
dimer. 

Altogether, the crystal structure of HsaDRS2 reveals 
that it shares the same global architecture with iicoDRS 
but differs by an enlarged catalytic groove, by a more 
electropositive surface potential and by a modified 
network of interactions in the dimer interface. 



Temperature stability of //saDRS2 versus £coDRS 

The sensitivity of both proteins to temperature has been 
compared using DLS and differential scanning fluorimetry 
(Figure 3 and Supplementary Figure S7). HsaDRS2 
unfolded at a temperature of >12°C below that of 
iicoDRS, in the 37°C^45°C range, under various biochem- 
ical conditions and in the presence of various additives. 
This result is supported by sets of synchrotron radiation 
circular dichroism spectra recorded at temperatures from 
20°C to 80°C (Supplementary Figure S8). For both 
proteins, a transition from an ordered to a disordered 
state occurs at temperatures similar to those observed in 
DLS and differential scanning fluorimetry. This low 
melting temperature was also valid for an HsaDRS2 
variant bearing different N- and C-terminal amino acids 
but not for HsaYRS2, which melted at ~55°C 
(Supplementary Figure S7). Significant stabilization was 
observed in the presence of Asp-AMS [(5'-0-[N-(L- 
aspartyl)sulfamoyl]adenosine; Glu-AMS, Gln-AMS, Tyr- 
AMS for other sulfamoyl-adenylates; AspOH, GluOH, 
TyrOH for amino alcohol-adenylates], a high affinity 
analog of the aspartylation reaction transition state (43), 
but neither in the presence of tRNA, adenosine triphos- 
phate or aspartic acid did (Figure 3 and Supplementary 
Figure S7). Asp-AMS had a specific effect at a 1:1 stoi- 
chiometry per monomer (not shown), suggesting that it is 
binding to the catalytic site and that this binding protects 
the protein from melting. Stabilization of aaRSs by 
binding of adenylate analogs has been taken advantage 
for successful crystallization (44). Herein, we observe sta- 
bilization against temperature induced melting. Along 
initial attempts to stabilize HsaDRS2 by site-directed mu- 
tagenesis, the replacement of amino acids of the hinge 
domain by those of EcoDRS led to a gain of 5°C in Tm 
(Supplementary Figure S7). 



Figure 1. Continued 

insertion and Arg271 (Arg222) in the motif 2 loop of class II 
synthetases. The 26° reorientation visualized by the movement of the 
axis of the external helix opens up the groove for the tRNA acceptor 
arm by 7-10 A in HsaDRS2 when compared with the tRNA bound and 
free states of iicoDRS. 
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Figure 2. Comparison of Z/.raDRS2 (left) and £coDRS (right) surface potentials. Surface potentials are depicted from positive to negative in blue to 
red with the same color scale. (A and B), tRNA-bound #.«?DRS2 (model) and EcoDRS (PDBid: 1C0A) with tRNA Asp (shown in purple). (C) The 
enzymes without tRNA. In (A), the complexes are observed along their two-fold symmetry axis. In (B and C), they are rotated as indicated to better 
visualize the interaction zone on the surface of the left monomer. The right monomer is colored with hydrophobic patches in green. 



Thermodynamics of tRNA binding to //s«DRS2 versus 
EcoDRS 

ITC (45) was used to gain insights into the thermo- 
dynamics of tRNA Asp binding to DRSs. In vitro tran- 
scribed tRNAs were used, as access to milligram 
amounts of native human mt-tRNAs remains out of 
scope (9). Mt-tRNA Asp distinguishes drastically from bac- 
terial tRNA Asp by primary, secondary and tertiary struc- 
tural elements (Supplementary Figure SI), leading to the 
absence of major aspartate identity elements and to a 
larger structural flexibility (25). Typical thermodynamic 
data are presented in Figure 4 and summarized in 
Table 3. Under the tested conditions, the affinities of 
HsaDRS2 for either mt- or Eco tRNA Asp were in the 
same range (K d ~250nM), whereas binding of EcoDRS 
to its cognate tRNA Asp was weaker by one order of mag- 
nitude (Kd of 3 uM) and binding to the non-cognate 
mt-tRNA Asp dropped the affinity by two orders of 
magnitude (K d of 22 uM). Enthalpy and entropy contribu- 
tions distinguished mitochondrial and bacterial cog- 
nate complex formation. 7/raDRS2/mt-tRNA Asp was 
accompanied by about twice as large contributions of 



both parameters than EcoDRS/Eco tRNA Asp formation. 
However, both values compensated each other to similar 
extents, leading to similar variations in free energy namely 
AG = —9.1 kcal/mol and — 7.5 kcal/mol, respectively. The 
difference between these free energies (A AG = — 1.6kcal/ 
mol) is in agreement with a AAG of —2.45 kcal mol -1 
measured by comparative kinetic properties of both 
aaRSs at the transition state (Table 1,3, and 4). 

Thermodynamic parameters for non-cognate complex 
formation (cross-complexes) revealed a pre-eminent role 
of the tRNA partner as compared with the synthetases 
(Table 4 and Figure 4A and B). Similarly, low binding 
enthalpies and entropies drive complexes involving the 
structurally stable Eco tRNA Asp regardless the enzyme. 
At opposite, binding of the unstable and highly flexible 
mt-tRNA Asp required far larger enthalpic and entropic 
compensations, regardless the enzyme. These large 
values are indicative of the adaptation of mt-tRNA Asp 
to both enzymes' surfaces involving likely conformational 
changes of the tRNA, releasing of water molecules along 
the binding surface and reciprocal adaptation of the 
nucleic acid and the protein. 
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Figure 3. Thermostability of HsaDRSl and EcoDRS. Melting tempera- 
tures as measured by DLS. (A) Thermostability of HsaDRS2 (dashed 
black line) compared with &oDRS (black line). (B) The presence of 
in vitro transcribed Hsa mt-tRNA Asp (light grey) or Eco tRNA Asp (dark 
grey) (in 1.2-fold molar excess) did not significantly influence the 
thermostability of either DRS. (C) Among the small substrates, only 
Asp-AMS (dashed grey line) influences the thermostability of the two 
enzymes positively [aspartic acid: black line, adenosine triphosphate 
(ATP): grey line, AspOH: light grey line]. Denaturation temperatures of 
duplicates performed with different enzyme preparations differed by <2°C. 



Thermodynamics of adenylate analog binding to 
HsaDRSl versus £coDRS 

Binding of the adenylate analog Asp-AMS to the two 
dimeric DRSs revealed two different binding sites within 
each DRS dimer (Figure 4 and Table 4). For both 
enzymes, a first adenylate bound to one monomer (n 
~0.5, i.e. 0.5 mol of adenylate per mol of monomer) with 
low affinity (Kd of 129 nM for HsaDRS2 and 29.5 nM for 
£coDRS). The second adenylate bound then to the second 
site with a higher affinity (Kd of 17.5 nM for HsaDRSl and 
of 2.95 nM for ZscoDRS), in favor of a positive-allosteric 
event. The gain in affinity for the second binding site was of 
similar amplitude (5- to 10-fold) in both enzymes, and 
affinities (in the nM range) were in agreement with the pre- 
viously measured K[ values for HsaDRSl (43). Binding of 
Asp-AMS to each monomer, however, followed a different 
thermodynamic route depending on the considered enzyme. 
For .EcoDRS, binding was both enthalpy- and entropy- 
driven. In contrast, binding to HsaDRSl was strongly 
enthalpy-driven. This different thermodynamic signature 
supports different binding mechanisms of Asp-AMS to 
HsaDRSl and to EcoDRS. 
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Figure 4. Isothermal titration calorimetry of HsaDRS2 and EcoDRS 
to tRNAs and to the adenyate analog Asp-AMS. (A) Titration of 
in vitro transcribed tRNAs by HsaDRS2. (B) Titration of in vitro 
transcribed tRNAs by &oDRS. (C) Titration of £foDRS and 
HsaDRS2 by Asp-AMS. 



DISCUSSION 

'Same-same but different' 

The crystal structure of HsaDRSl revealed a typical 
overall bacterial-type DRS architecture (30, 38-40) that 



superimposes with the structure of EcoDRS with a 
global rmsd <2 A (Figure 1). The two functionally import- 
ant domains, catalytic site and anticodon-binding domain 
overlap with an rmsd of < 1 .5 A. The important differences 
between the two structures are (i) the angle formed by the 
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Table 3. Comparative thermodynamic properties of //.«/DRS2 and £coDRS: thermodynamic parameters of DRS interaction with tRNA, as 
measured by ITC 

DRS + tRNA(transcripts) n a K A (uM) AH (kcal moP 1 ) -TAS (kcal moP 1 ) AG (kcal mor 1 ) 

HsaD~RS2 + Hsa mt-tRNA Asp 0.48 ± 0.26 0.26 ± 0.18 -20.3 ± 0.7 +11.2 ±0.4 -9.1 ±0.5 

HmDRS2 + Eco tRNA Asp 0.87 ± 0.16 0.24 ± 0.01 -14.0 ± 0.7 +5.0 ± 0.7 -9.0 ± 0.5 

iicoDRS + Eco tRNA Asp 0.76 + 0.01 3.1 ± 0.1 -13.0 ± 3.1 +5.5 + 3.2 -7.5 + 0.5 

EcoDRS + Hsa mt-tRNA Aspb 0.39 + 0.02 22.0 -30.8 +24.4 -6.4 

All data were obtained at 25°C (T = 298°K) and correspond to mean values from at least duplicated experiments with independent and freshly 
prepared enzyme and tRNA lots. 

" l n stands for the number of moles of Asp-AMS bound per mole of DRS monomer. 
b Cross-binding between &oDRS and H.sa mt-tRNA Asp performed only once. 



Table 4. Comparative thermodynamic properties of Z/.raDRS2 and £coDRS: thermodynamic parameters of DRS interaction with Asp-AMS, as 
measured by ITC 

DRS + Asp-AMS n a K A (nM) AH (kcal moP 1 ) -TAS (kcal MoP 1 ) AG (kcal moP 1 ) 

i/.?flDRS2-monomer 1 0.48 + 0.02 129 + 38.0 -13.18 + 0.39 +3.9 + 0.78 -9.38 + 0.18 

//.s«DRS2-monomer 2 0.29 + 0.04 17.5 + 12.0 -21.8 + 0.31 +10.57 + 0.71 -10.63 + 0.44 

£coDRS-monomer 1 0.49 + 0.02 29.5 + 5.0 -5.50 + 0.01 -4.75 + 0.13 -10.24 + 0.10 

EcoDRS-monomer 2 0.44 + 0.01 2.95 + 1.34 -8.16 + 0.35 -3.52 + 0.61 -11.63 + 0.28 

All data were obtained at 25°C (T = 298°K) and correspond to mean values from at least duplicated experiments with independent and freshly 
prepared enzyme and tRNA lots. 

' l n stands for the number of moles of Asp-AMS bound per mole of DRS monomer. 



prokaryal-specific insertion domain and the catalytic 
domain, which is more opened by 26° in the mt-aaRS; 
(ii) the overall electropositive surface potential (as 
observed for other mammalian mt-aaRSs); and (iii) the 
dimer interface, which is weaker in terms of salt-bridges 
and hydrogen bonds. 

A distinguishing property between the two proteins is the 
unprecedented observation for a human mt-aaRS to melt 
at near body temperature (37°C) and ~12°C lower than 
iicoDRS, in line with a restricted thermo-tolerance of the 
translation and replication processes in mitochondria (46). 
It is likely that the weakened dimer interface is one of the 
contributors to this fragility. Fine prediction of the molecu- 
lar reasons for the melting temperature of a protein is 
difficult and may be dependent on as little as a single 
amino acid (47). Interestingly, a thermosensitive mutant 
of EcoDRS in which proline 555 was affected, has been 
reported (48). The equivalent position in HsaDRS2 has 
kept an equivalent proline residue (Pro 605), so that the 
contribution of this position to thermosensitivity is ruled 
out. We favor the contribution of the hinge domain, as 
partial stabilization by mutation in HsaDRS2 was 
successful. 

ITC was used here for the first time on tRNA/synthetase 
couples to examine thermodynamics of cognate 
mt-tRNA Asp and non-cognate Eco tRNA Asp interactions. 
When tested under strictly the same conditions, HsaDRS2 
binds both its cognate and non-cognate tRNA sp with 
a same sub-micromolar Kd, whereas EcoDRS binds 
cognate Eco tRNA Asp in the low micromolar range and 
non-cognate mt-tRNA Asp in the high micromolar range. 
It is anticipated that the parameters would be closer to, 
but still distinct, from those for the bacterial cognate 



pair when measured with native post-transcriptionally 
modified tRNAs. Indeed, modifications are known to 
stabilize structures (49), but the modification pattern of 
mt-tRNA Asp is far beneath that of Eco tRNA Asp (see 
Supplementary Figure SI). Binding between cognate mito- 
chondrial partners requires large enthalpic and entropic 
compensations at opposite to binding between the 
cognate bacterial partners where compensations were 
small. Cross-binding revealed that this difference is 
mostly linked to the tRNA. Binding of the floppy 
mt-tRNA involved large positive enthalpic and large 
negative entropic contributions, whereas binding of the 
well-folded and stable Eco tRNA involved low enthalpy/ 
entropy compensations. These results are in line with 
thermodynamic parameters determined for cognate and 
non-congate tRNA/synthetase couples measured earlier 
by fluorescence titrations and fast kinetic techniques 
(50,51). 

Altogether, comparative investigation of biophysical 
properties of HsaDRS2 and £coDRS shed new light on 
the long-term open question as to non-reciprocal amino- 
acylation of bacterial and mitochondrial tRNA/ 
synthetase pairs (24). They revealed that the human 
enzyme has gained structural adaptability to tRNA sub- 
strates by combining several properties. A more electro- 
positive tRNA-binding surface [see also (52)] and an 
opening of the cleft between the insertion domain and 
the catalytic domain by an angle of as much as 26° 
improve accessibility and binding of the tRNA end to 
the catalytic site. Strong binding likely reflects an evolu- 
tionary adaptation of the mt-aaRS to compensate for 
high flexibility of the mt-tRNA and for the absence of 
otherwise important identity elements. Whether the 
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Figure 5. Molecular surface potential of bacterial (E. coli or Thermus 
thermophilus) and mammalian (bovine or human) mitochondrial aaRSs. 
On the left, bacterial complexes (PDBids: 1C0A for £coDRS/tRNA, 
1SRS for 77/;SRS/tRNA, 2IY5 for 77/iFRS/tRNA and 1H3E for 
77/)YRS/tRNA) are shown on the left with monomers A in blue, 
monomers B in green and tRNAs in pink, indicating the binding site 
of the cognate substrates. In the case of tetrameric bacterial FRS, 
monomers C and D are depicted in cyan. The electrostatic surface 
potential of free aaRS forms is represented in the same orientation as 
the complexes for bacterial (1EQR for £coDRS, 1SRY for 77/iSRS, 
1B7Y for 77/iFRS, 1H3F for TthYRS without its C-terminal domain) 
and mitochondrial enzymes (this study for i/.s«DRS2, 1WLE for 
BfoSRS2, 3TUP for HsaFRSl, 2PID for HraYRS2). Overall, the 
tRNA binding surface seems to be more electropositive in mitochon- 
drial forms, even for FRS, which has a different structural organization 
(monomer instead of heterotetramer). 



binding step is accompanied by a concerted movement 
of the insertion domain towards the catalytic site to 
further lock the tRNA, is unknown. The larger structural 
plasticity of the mt-aaRS (enabling it to deal with 
the highly flexible and degener ated mt-tRNA Asp ) does 
not harm its interaction with classical and less flexible 
Eco tRNA Asp . At opposite, the bacterial aaRS is 
unable to bind properly the mt-tRNA. Binding to its 
cognate more rigid tRNA sp is closer to a lock-and-key 
type process. 

Gain of plasticity as reflected here by unusual thermo- 
dynamic properties, may be a more general property of 
mt-aaRSs that have to deal with degenerated mt-tRNAs. 
An extreme case of plasticity has recently been described 
with the human mitochondrial FRS that undergoes 



dramatic conformation change on tRNA binding (21). 
The near-body thermal unstability of HsaDRS2 as well 
as its electropositive surface potential suggest that differ- 
ent parameters may stabilize the protein in vivo. These 
may be of biochemical nature or correspond to co-factors, 
partner proteins or peculiar sub-mitochondrial localiza- 
tion, as is the case for mt-ribosomes and mt-elongation 
factor that bind to the inner membrane (9). 

Cooperative allostery in class II dimeric DRSs 

A long-standing question for dimeric aaRSs concerns the 
communication mechanisms between the two catalytic 
sites (13). Data obtained herein by ITC demonstrate un- 
ambiguously that both iicoDRS and HsaDRS2 undergo 
an allosteric event on binding to the adenylate analog. 
The first monomer binds the analog with a weaker 
affinity than the second monomer. Thermodynamic signa- 
tures strongly suggest a different binding mechanism for 
each enzyme: interaction with HsaDRSl is a strongly 
exothermic process, characterized by a large negative 
enthalpic term, whereas binding to EcoDRS is associated 
with moderate enthalpic and entropic contributions. 
These differences likely reflect large obligatory rearrange- 
ments in HsaDRS2 but not in £coDRS. 

Allostery has been reported for several aaRSs (53-56). 
Half of the site activity was mainly observed. Herein, we 
report direct cooperative allostery for small substrate 
binding, as detected by ITC. Previous ITC investigations 
performed on other aaRSs had not detected cooperative 
allostery (53-56). 



CONCLUSION AND PERSPECTIVES 

Resolution of the crystal structure and exploration of bio- 
physical properties of a human mitochondrial aaRS, 
revealed similarities but also distinctive features from a 
prokaryotic homolog: a conserved architecture and a 
conserved cooperative allosteric communication within 
the dimers but also increased structural plasticity. We 
propose that the modification of intrinsic biophysical 
properties and new surface properties, along with conser- 
vation of size and structural framework, represents an 
evolutionary mechanism for adaptation of nuclear- 
encoded proteins to degenerated mitochondrial-encoded 
RNAs. This mechanism is alternative to domain extension 
and sustains adaptation of two genomes subjected to dif- 
ferent evolutionary constraints. The possibility that 
mt-aaRSs also undergo functional expansion by structural 
inventions and have partner proteins as described for 
cytosolic aaRSs (57,58) is not excluded but remains to 
be tackled. The properties of an mt-aaRS reported 
herein provide new insights towards the understanding 
of the growing number of human disorders linked to mu- 
tations in the corresponding genes (9,59,60). Deciphering 
the fine molecular triggers leading to both specific (i.e. 
plasticity) and shared (i.e. allostery) properties of human 
mt-aaRSs, and their evolutionary related homologs 
requires delineation of discriminating and 
non-discriminating amino acids within the aaRS 
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architecture. Bioinformatic tools enabling such an analysis 
are under development. 
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